Renormalization group approach to nonextensive statistical mechanics 
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We analyze a simple classical Hamiltonian system within 
the hypothesis of renormalizability and isotropy that essen- 
tially led Maxwell to his ubiquitous Gaussian distribution 
of velocities. We show that the equilibrium-like power-law 
energy distribution emerging within nonextensive statistical 
mechanics satisfies these hypothesis, in spite of not being 
factorizable. A physically satisfactory renormalization group 
emerges in the {q,T q ) space, where q and T q respectively are 
the entropic index characterizing nonextensivity, and an ap- 
propriate temperature. This scenario enables the conjectural 
formulation of the one to be expected for d-dimensional sys- 
tems involving long-range interactions (e.g., a classical two- 
body potential oc r~ a with < a/d < 1). As a corollary, we 
recover a quite general expression for the classical principle of 
equipartition of energy for arbitrary q. 

PACS number(s): 05.10.Cc; 05.20.-y; 05.20.Gg; 05.70.Fh 



If we allow ourselves to use a contemporary vocabu- 
lary, we may say that the essential hypothesis that led 
James Clerk Maxwell to the classical distribution of ve- 
locities that is named after him were two, namely renor- 
malizability and isotropy. Let us be more specific. If wc 
have the one- and two-degrees-of-freedom (hamiltonian- 
like) quantities h'(x;a') = a'\x\ z and h(x,y;a,b) — 
a\x\ z + b\y\ z , is there any function f(h) such that ex- 
act renormalization occurs upon reducing the degree of 
freedom yl More explicitly, is there any f(h) such that 
dy f(h) oc f(h')l It is well known that Maxwell 
thought of the exponential form, i.e., f(h) — exp(-h). 
Isotropy of course fixed z — 2 (consistently with New- 
tonian mechanics). The fact that factorization occurs 
(i.e., f(h(x,y;a,b)) = f(h'(x;a))f(h'(y;b))) can be con- 
sidered, from the present viewpoint, as a simplifying 
aside consequence. Incidentally, it is worthy to note that, 
at Maxwell's time, his arguments constituted a break- 
through; indeed, the distribution of velocities occasion- 
ally employed at that time was flat within some interval, 
and zero outside. 

At this stage, let us point out that the exponential 
function is but the q — 1 member of an entire fam- 
ily of functions, namely the g-exponentials |l|-|3| f(h) = 
exp q (—h) = [1 — (1 — q)h]~ which, for a continuous 
range of values of q, are renormalizable, and which also 
can satisfy isotropy. Let us more precisely define the q- 
exponential function: 



[i + (i- q )t]— (ten-, qeTZ) 



(1) 



For q < 1, the q-exponential function vanishes for t < 
— 1/(1 — q) and continuous and monotonically increases 
from to oo when t increases from —1/(1 — q) to oo. For 
q > 1, the ^-exponential function continuous and mono- 
tonically increases from to oo when t increases from 
— oo to l/(q — 1), remaining divergent for t > l/(q — 1). 
We easily verify that e* = lim g ^i + o e q = lim g ^i_o e* = 
e* (Vt). 

Let us now illustrate, for a simple case, the above men- 
tioned renormalizability: 



dx2 exp q (— ai\xi\ zi — a 2 \x 2 \ Z2 ) 

) 

= A 2 exp q ,(-a' l \x 1 \ Zl ) (ai, a 2 , z x , z 2 > 0) (2) 



where 
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l-q> 1-q 
a' 1 (l-q')=a 1 (l-q) 



(3) 
(4) 



and 



A 2 = 



d£ exp q {-a 2 \£\ z >) = 

(i- g )V, 2 r(i+i + T ^) 
1 

V»-l '2 ) 



for q < 1 
for 3=1 
for q > 1 



(5) 



This type of renormalization remains exact for the follow- 
ing A-degrees-of-freedom dimensionless Hamiltonian: 



JV 



h N {{xi}; {a,}) = y^^ai\xi\ Zi 
i=l 

We straightforwardly verify that 

dx N exp q (-h N ({x. t }; {a,})) 
= A N exp <3 /(-/ijv-i({a: J }; {ai})) 



(6) 



(7) 



where q' and An are respectively given by Eqs. (3) and 
(5) , with z/v and a at replacing z 2 and a 2 . Eq. (4) is gen- 
eralized into a-(l — q') = aj(l — q) (i = 1, 2, N — 1). 
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It is clear that there is no major difficulty in reducing, at 
every step, more than one degree of freedom. Hamilto- 
nian (6) with Zi = 2 (\/i) corresponds, of course, to N free 
particles in one dimension, or N / d particles in d dimen- 
sions; isotropy is automatically satisfied in such cases. 
Also, with Zi — 2 for N/2 degrees of freedom and Zi — z 
for the other N/2 degrees of freedom, it corresponds to 
N/2 one-dimensional anharmonic oscillators. 

Let us now make the junction of the above renormal- 
izability ideas with the nonextensive statistical mechan- 
ics introduced a decade ago further implemented in 
||U and applied in a considerable variety of systems, 
such as Levy and correlated H anomalous diffusions, 
peculiar velocities in spiral galaxies 0], turbulence in 
electron plasma j[o|, fully developed turbulence H), ci- 
tations of scientific papers fl^| , reassociation in folded 
proteins G3], quantum entanglement fllfl , and others 
(for recent reviews, see ||[L5|]). To be more precise, 
this formalism addresses systems which, in one way or 
another, exhibit spatial and/or temporal long-range in- 
teractions, multifractal structures, dissipation, and re- 
lated anomalies. The basic formal standpoint for thermal 
equilibrium (or equilibrium-like stationary states) con- 
sists in optimizing, under appropriate constraints, the 
entropic form S g = [1 - J2iPi}/il ~ 1] (Q e ^ {Pi\ is 
the set of probabilities associated with the microscopic 
states). In particular, for the canonical ensemble, the 
constraints are Q J2iPi = 1 an d J2iPiEi = U q ({Et} is 
the set of the energies of the microscopic states), where 
Pi = Pi/J2jPj- ^ i s cl ear that the q — s- 1 limit repro- 
duces the standard Boltzmann-Gibbs recipe; in particu- 
lar Si — — J2iPi^ n Pi ( we are using fcs = 1). Also, if A 
and B are two probabilistically independent systems, we 
have that S q {A+B) = S q (A)+S q (B) + (l-q)S q (A)S q (B), 
hence q = 1, > 1, < 1 respectively correspond to exten- 
sive, subextensive and superextensive systems. 

Optimization for arbitrary q straightforwardly leads || 

to 



exp q (~(Ej - Ug)/T q ) 
Z„ 



(8) 



where Z q = £\ exp q (-{E j -U q )/T q ) with T q = Tj^jPj, 
1/T = (3 being the Lagrange parameter associated with 
the energy constraint. As we see, the similarities with 
Boltzmann-Gibbs statistical mechanics are quite strik- 
ing. For instance, it can be verified that 1/T = dS q /dU q . 
Also, systematically, the function e q and its inverse 
ln q t = (i 1-9 — 1)/(1 — q) play the usual roles of the 
exponential and logarithm ones. In particular, F q = 



U q - TS q 



-Thi q Z q and U q = — <9(ln g Z q )/d(3, where 



Z q is defined through ln g Z q = ln 9 Z q — j3U q . 

The application of this formalism to the classical 
Hamiltonian (6) yields 



Z„ 



d x exp q {-{h N ({xi\; {a,}) - U q )/T q ) 



I q>{ ^ } eXPq(Uq/T q ), 



(9) 



[l + (l-q)U q /T q ]T q and I q , {bi} 



where r q \ 

J d N x exp q (—h N ({xi};{bi})) is easily, though tediously, 
calculated from the recursive application of Eqs. (2-5). 
Analogously we obtain 



Z q = I d N x [expq{-(h N {{x l }-{a i \)-U q )/T q )} q 



I qd s^ } [exp q (U q /T q )] q , 



(10) 



Zn = Z a we 



where q = 2 — l/q. By using the identity 
can solve the present set of implicit equations. It follows 
the generalized classical equipartition principle (see also 
@) 



U a 



N 1 



To 



(11) 



Consequently r q 



T q . It is kind 



i + (i-g)£i=i 

of remarkable that such a strongly Boltzmann-Gibbs-like 
equality as Eq. (11) does hold, in spite of the fact that 
the involved distributions are not factorizable. 

Finally we must focus on the probability distributions: 

P^ N) oc (p q N) y (x exp s (-^f-\ Then, by reducing 

one degree of freedom (as done in Eqs. (2-5)), we obtain 



Pf-V - J dx N PW 

oc exp (2 -i/ q >) (~q h N -x ({x,} ; {a l /r q })) . (12) 
Consequently, by considering z% = z (Vi) we obtain 
, 1 - q + qz 



and 



T' 

g 



1 — q + z 



1 - q + qz 



(13) 



(14) 



where we have used the fact that T q oc r q with a propor- 
tionality coefficient which preserves, as easily verifible, 
the initial value of q, and not the renormalized one. The 
results corresponding to say N anharmonic oscillators are 
completely analogous. 

The recurrence for q has only one (double) fixed point, 
namely q — 1. This point locally is an inflexion one, but 
globally is an attractor (like the tangent bifurcation of 
the logistic map at the entrance of the cycle-3 window): 
see Fig. 1. In other words, Boltzmann-Gibbs statisti- 
cal mechanics is, as well expected, the correct one for 
the independent-particle system we have focused on here. 
We also verify that (i) T q = is an invariant subspace of 
the renormalization group; (ii) there is no T g -flow at the 
q = 1 subspace, and (iii) Eq. (14) is invariant through 
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the transformation T q — > AT g (A > 0) (which means that, 
without loss of generality, we can always start the recur- 
rence with say T q = 1). A typical two-dimensional flow 
is shown in Fig. 2. Before going on, let us mention that, 
if we start the recurrence with < q < 1, the flow in 
the q, T q space is smooth and monotonic (q approaches 
unity, and T q approaches zero). If we start with q not 
much larger than unity, the flow is smooth and mono- 
tonic (both q and T q increase) as long as the iterations 
provide q-images not exceeding 1 + z (divergence of q' 
in Eq. (13)). After that, the behavior becomes phys- 
ically meaningless. Indeed, the g-image is quickly sent 
to quite negative values, then approaching q = 1 from 
below. Concomitantly, T q changes sign, and approaches 
zero from below. This mathematical artifact is easy to 
understand; indeed, the integrals involved in the renor- 
malization are not defined for q > 1 + z. 

Since the renormalization takes N into N — 1, all fixed 
points necessarily correspond to the thermodynamic limit 
N — > oo. The same type of flow is expected to correspond 
to any (classical or quantum) system with noninteracting 
elements. If there are collective short-range interactions 
(e.g., a classical <i-dimensional fluid with two-body inter- 
actions decaying as r~ a , or an Ising ferromagnet with 
coupling constant decaying as r~ a , with a > d), we still 
expect q — 1 to be the unique (globally attractive) fixed 
point, but on the q = 1 invariant subspace there might 
be a finite T q unstable fixed point, T q — and T q — » oo 
being stable fixed points, respectively corresponding to 
the ordered and disordered phases. However, if the in- 
teractions are long-range ones, we expect the double fixed 
point q = 1 to split into two fixed points, namely a glob- 
ally attractive one at q = 1, and a globally repulsive one 
at q* > 1, whose value should depend on (a, d) (natu- 
rally, q* is expected to continuously approach unity when 
a/d — > I). See Fig. 3. A clarification is needed: why we 
rather expect q* > 1 and not q* < 1? This comes from an 
everyday increasing evidence in a variety of noncxtensive 
systems (electron-positron annihilation Jl7| , quark-gluon 
plasma p8|| , granular matter [[ljj]], cosmology d = 1 
system of inertial classical planar rotators ferromagncti- 
cally coupled at long distances among others), where 
the Boltzmann exponential distribution of energies is re- 
placed by a long-tailed power-law, and this precisely is 
what occurs for q > 1 (see Eq. (8)). One more clarifi- 
cation is needed: why we rather expect the q = 1 fixed 
point, and not the q ^ 1 one, to be the (globally) at- 
tractive one? Once again, this comes from increasing ev- 
idence 22 1 that the Boltzmann regime is the ultimately 
stable one for any finite-size system (it corresponds to 
the ZzmAr^oo/jmt^oo), whereas the nonextensive regime 
(i.e., q ^ 1) emerges in the limt—>oolimN—>oo ordering (see 
also Fig. 4 of Jl5| ). Moreover, the present conjecture is 
consistent with the topology recently found in |^3| for 
anomalous diffusion, namely a stable fixed point at q = 1 
and an unstable one at q = 2. 



Summarizing, we have shown that the renormalizabil- 
ity and isotropy principles which guided Maxwell argu- 
ments in establishing his celebrated Gaussian distribu- 
tion of velocities is not a privilege of the exponential 
function, but it is shared by an entire power-law family of 
functions (the g-exponentials), which includes the stan- 
dard exponential as the q = 1 limiting case. This pro- 
vides a remarkable bridge with nonextensive statistical 
mechanics, where precisely the g-exponentials play a cen- 
tral role (just as the exponential does within Boltzmann- 
Gibbs statistical mechanics). This observation straight- 
forwardly enabled the establishment of the g-generalized 
classical principle of equipartition of energies. It also en- 
abled the formulation of an exact renormalization group 
in the (q, temperature)- space, which shows that the 
globally stable fixed point q = 1 (i.e., Boltzmann-Gibbs 
statistical mechanics) is a double one. This allowed the 
conjecture that, for systems including long-range interac- 
tions, this fixed point splits into two, one of them (q = 1) 
remaining globally stable, and the other one (presumably 
q > I) being globally unstable. This situation respec- 
tively reflects the fact that the extensive or nonextensive 
thermostatistical are to be observed in the limits where 
the first to achieve infinity is the time or the size. 

One of us (RSM) acknowledges S. Matteson and P. 
Grigolini for warm hospitality at the Department of 
Physics/UNT. 
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FIG. 1. Noninteracting Hamiltonian (with z — 2): Renor- 
malization of q (any initial value of q is ultimately attracted 
by q = 1; in particular, if we start with q > 1, the successive 
values of q might increase for a few steps, but then the im- 
age is sent to the region q « 1, and then it monotonically 
increases up to q = 1). Renormalization alternates between 
the curves q' = q and q' given by Eq. (13), which exhibits a 
vertical asymptote at q = z + 1. 



FIG. 2. Noninteracting Hamiltonian (with z = 2): Renor- 
malization group flow in the {q,T q ) space, by starting with 
three different initial values of q; the initial value for T q is 
always taken to be unity (see the text). 



FIG. 3. Typical long-range interacting Hamiltonian: Con- 
jectural projection on the q-axis of a two- or more-dimensional 
renormalization group flow. Two fixed points are expected, 
namely at q = 1 (attractor) and at q — q* (repulsor); q* 
should approach unity when the range of the interactions de- 
creases (for all a above some critical value, q* is expected to 
remain equal to unity). 
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